Predicting plastic flow events in athermal shear-strained amorphous solids 
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We propose a method to predict the value of the external strain where a generic amorphous solid 
will fail by a plastic response (i.e. an irreversible deformation), solely on the basis of measurements 
of the nonlinear elastic moduli. While usually considered fundamentally different, with the elastic 
properties describing reversible phenomena and plastic failure epitomizing irreversible behavior, we 
show that the knowledge of some nonlinear elastic moduli is enough to predict where plasticity sets 
in. 



Introduction: Studies of plasticity in amorphous 
solids have always been hampered by the lack of any 
method to identify a-priori the locations of dissipativc 
(plastic) events. This is in sharp contrast with crystals, 
where plasticity can be assigned to the motion of iden- 
tifiable, topological defects like dislocations, the discov- 
ery of which has triggered tremendous theoretical break- 
throughs. In the absence of any clear-cut definition of 
"shear transformation zones" , studies of amorphous sys- 
tems are thus still in search for reliable predictors of yield- 
ing. Obvious observables such as local stress or density 
have proven unreliable [l|, [If. More promising are studies 
of local elastic fluctuations ■ It appeared that plastic 
failure correlates with soft elastic regions. This could be 
expected as yielding involves the crossing of saddle points 
at which the shear modulus vanishes before it presents a 
singular behavior Q. Nevertheless the values of the lin- 
ear elastic moduli alone do not carry predictive power; 
we cannot just say that failure will occur in any given 
softer region in space. The search for reliable predictors 
of plasticity has thus become a major issue in studies 
of amorphous solids. In this Letter we will show that 
an accurate predictor of plastic failure in an athermal 
amorphous solid can be constructed, as soon as higher 
order derivatives of the potential function are involved. 
Our findings not only offer a predictive tool for the onset 
of failure, but also point out the importance of nonlin- 
earities, and in particular those that couple nonlinearly 
"softening" regions with strain at larger scales. 

To fix ideas, imagine a simple shear deformation ap- 
plied to a given piece of amorphous solid (for simplicity 
in 2D, with immediate extensions to 3D). A small strain 
increment Sj corresponds to a change of the i'th parti- 
cle positions — > r[ as: x[ — Xi + Sjyi, y[ = yi , In 
athermal quasi-static conditions (T — > 0, 7 — > 0), the 
system lives in local minima, and follows strain-induced 
changes of the potential energy surface @, Therefore, 
the particles do not follow homogeneously the macro- 
scopic strain, and their positions change as n — s- r ■ + 
where Ui denotes non-affine displacements. Around some 
stable reference state at 7 = 70, the field «i, the system 
energy, and internal stress a xy are smooth functions of 
7. We can thus write: 
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FIG. 1: A typical stress vs. strain curve in a system of 4096 
particles in two dimensions obtained in the athermal limit. 
Every elastic (reversible) increase in stress if followed by a 
sudden plastic (irreversible) drop in stress. The aim of this 
Letter is to predict the value of the strain at which the next 
plastic drop will take place. 

As the strain increases, the system must eventually lose 
mechanical stability; the "elastic branch" on the stress 
curve ends in a discontinuity as the system fails via a first 
subsequent "plastic event" , see Fig. Q] It is precisely at 
this instability, say at 7 = 7p, that the function o~ X y(l) 
loses its analyticity. Accordingly we recognize that the 
radius of convergence of the series (|TJ) is precisely \"fp — 
7o|, where jp can be larger or smaller than 70. 

Our method to predict where plasticity sets in rests on 
two ideas. The first is that the coefficients in Eq. ([T]) 
contain relevant information about the value of 7 where 
analyticity is lost. This comment is reminiscent of high 
temperature expansions in critical phenomena, where the 
knowledge of a substantial number of derivatives (and 
customarily using some Pade resummation) can shed im- 
portant light on the nature of the critical point [Toj • The 
second is that, although we do not have access to a large 
number of derivatives (see below), we actually know [7| 
the nature of the singularity at 7p, and we can use this 
knowledge to reach predictability which otherwise would 
be out of hand. Recall that as long as the system remains 
in mechanical equilibrium (i.e. along the elastic branch) 
the force fi on every particle are zero before and after an 
infinitesimal deformation; In other words @, [H| with U 
the potential energy 
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This condition introduces the all-important Hessian ma- 
trix ffy and the 'non-affinc force' Sj which can both be 
computed from the interparticle interactions. We rewrite 
this condition as 
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where the second equation results from expanding in the 



eigenfunctions of H, Hijipj = X^f] the last esti- 
mate stems from our knowledge that in finite systems the 
plastic event is associated with a single eigenvalue going 
through zero when the systems slides over a saddle. De- 
note the critical eigenvalue as Ap. Eq. (Q| integrates to 
provide the distance of the non-affine field Ui from its 
value at j P , Ui(j) - Ui(jp) = X(7)i/>,| P) , where X(j) is 
a function of 7 only, satisfying X(j P ) — and 
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Finally, we use the crucial assumption [?J that the eigen- 
value Ap crosses zero with a finite slope in the X- 
coordinate system itself, where distances are measured 
along the unstable direction: 



Ap pa AX + 0{X 2 
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Together with Eq. ([5]) and asserting that Bj is not sin- 
gular (it is a combination of derivatives of the potential 
function [HI]), implies that 
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These results are now used to determine the singularity 
of the stress at jp . We start with the exact result for the 
shear modulus 0, HH 



A 4 = Ms — B ■ H.~ x ■ B/V , 
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where /ip is the Born term. Using Eqs.© and we 
conclude that near 7p we can write the shear modulus 
as a sum of a regular and a singular term, 
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We thus assert that in the vicinity of 7p: 
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, 7 < 7p . 
(10) 

Indeed, the stress begins to go down before the plastic 
even takes place, but this is not seen in Fig. [T] since this 
happens very sharply. We can now come back to the 
question of predicting plasticity by looking at derivatives 



at any given point 70 < 7p. Near 70, we know that 
&xy ~ m(7 ~ To) + ■ • ■ • Using pU|) we write the ansatz 

Oxyh) = ffo + M(7-7o) + aVTP _r 7 + &(7P-7) 3/2 j (H) 
from which we can recalculate the derivatives at 7 = 70: 
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Note that the shear modulus has disappeared from these 
expressions, becoming irrelevant for the subsequent pre- 
diction of the instability threshold. By measuring these 
three derivatives at 70 we can determine all the unknowns 
in Eq. (jlll) . 7p solves a 4th order polynomial, and 
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is the only physical solution among the four available 
ones. Next we test these results in a specific model. 

Model and numerical procedures: Below we em- 
ploy a model system with point particles of two 'sizes' 
but of equal mass m in two-dimensions, interacting via a 
pairwisc potential of the form 
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where ry is the distance between particle i and j, £ is the 
energy scale, and x c is the dimensionless length for which 
the potential will vanish continuously up to q derivatives. 
The interaction lengthscale Xij between any two particles 
i and j is Ay = 1.0A, Ay = 1.18A and Ay = 1.4A for two 
'small' particles, one 'large' and one 'small' particle and 
two 'large' particle respectively. The coefficients C21 are 
given by 
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We chose the parameters x c — 7/4, k = 10 and q = 6. 
The unit of length A is set to be the interaction length 
scale of two small particles, and e is the unit of energy. 
The density for all systems is set to be N/V = 0.85A -2 . 
We employ an athermal quasi-static scheme which con- 
sists of imposing an affine transformation to each particle 
of a configuration, followed by a potential energy mini- 
mization under Lees-Edwards boundary conditions [l3| . 
In this scheme one can obtain purely elastic trajectories 
of stress vs strain 12j, which allows for the calculation 



of the total derivatives of stress with respect to strain 
using finite differences; we choose the stopping criterion 
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FIG. 2: Color online. Upper panel: Straining from an elas- 
tic branch of the elasto-plastic steady state in a system of 
N = 484 particles. Shown are the stress vs. strain curve up 
to the first plastic event (blue rhombi) and the value of "/p 
as predicted by Eq. (|14p for every value of 7 (green squares) . 
In addition we show the prediction of the simplistic Taylor 
expansion (|17[1 (red circles). Lower panel: the same measure- 
ments starting from the equilibrium isotropic stress-less state. 
Note the region where the even derivatives build up. 

for the minimizations to be |Vj£/| < 10~ 25 -| for every 
coordinate Xi , and select the strain increment for taking 
derivatives to be £7 = 5 x 10~ 7 . 

It is useful to compare our prediction with a simple 
failure criterion which uses a naive Taylor series of the 
form cr xy (j) = (To + ^(7 - 70) + ^B 2 (j - 7o) 2 + 5-83(7 - 
7o) 3 + ^484(7 — 7o) 4 - From here 7p is estimated as the 
point at which |^| 7=7P = 0-the difference between this 
point and the actual yield point is negligible compared 
to all other strain scales entering the problem This 
is the solution of the equation: 

M + S 2 (7p-7o) + 5-B 3 (7p-7o) 2 -l-i J B4(7p-7o) 3 =0 . 

(17) 

Predictions: Fig. [5] in the upper panel demonstrates 
the quality of our predictor for a typical initial condition 
taken in the elasto-plastic steady state in a system with 
N = 484. Here, 70 corresponds to the strain at which the 
elastic branch was reached via a prior plastic event. The 
total length of the branch on the stress-strain curve is 
thus exactly 7p — 70- We show the stress vs. strain curve 
in (symbols) and the values of 7p predicted using cither 
Eq. (fl4"|) or (IT71) at each value of the external strain. We 
see that the naive Taylor expansion fails except right at 
7p, whereas our accelerated method works very well over 
nearly the whole length of the elastic branch. Of course, 
in the elasto-plastic steady state the strain range between 
plastic events is rather limited-as reflected by the scale of 



abscissa. But our data shows that the transition between 
one plastic event and the next in entirely controlled by 
a single mode, i.e. by a single soft zone. We will shortly 
come back on this issue. 

We also test our predictor starting from an equilibrium 
isotropic state obtained by thermal annealing. Here we 
need to predict the first plastic event which occurs upon 
increasing the external strain. This situation could be ex- 
pected to be more tricky since at 7 = 70 = all the even 
derivatives d 2n a xy /d / y 2n vanish on the average. We could 
expect that, upon straining, these derivatives must build 
up before predictability of plasticity is achieved. In fact, 
as illustrated on the lower panel of Fig. [5] predictability 
is often achieved rapidly, even though the strain range is 
considerably larger than in the upper panel. In steady 
state the vast majority of the clastic branches end up in 
a plastic event whose onset can be predicted as demon- 
strated here. In contradistinction, the elastic branch em- 
anating from the isotropic equilibrium state may exhibit 
large non-affine elastic events which resemble a typical 
precursor to a plastic failure but avoids it by eventu- 
ally stabilizing. Our derivatives will pick up these elastic 
events and will incorrectly predict a plastic failure before 
landing on the right prediction. This scenario is shown 
in Fig. [3] in which two non-affme large events occur be- 
fore the first plastic event. All three derivatives undergo 
huge changes at the non-afiine events and then dive to 
—00 at the true plastic event. This seems to be a trace 
of the typical presence of several interacting eigenmodes, 
not necessarily well-aligned with the direction of shear- 
ing. We stress that this phenomenon is quite rare and 
the normal situation is the one that is exhibited in Fig. 
[2] lower panel. 

Dependence on system size: The data shown in 
Figs. [2] and [3] pertain to a relatively small system of 484 
particles. We therefore must raise the important question 
of how the predictability of the plastic failure depends 
on the system size. To study the size dependence of the 
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FIG. 4: Log-log plot of the range of strain |7 — 7p| where 
S'yp/Ayp < 0.15 as a function of the system size in the elasto- 
plastic steady state. The slope of the continuous fit line is 
77 = —0.61. The inset shows how the measurement is done, 
< 0.15. 

predictability we measured the distance Ajp = |7 — 7p| 
for which the error in estimating 7p, denoted as 5"/p, 
satisfies 5jp/A~fp < 0.15, for systems of varying sizes 
(see inset in Fig. 0}. The results are shown in Fig. 0] 
indicating that the range of predictability reduces like a 
power law with the system size, 

A 7P ~jsr> , t) w -0.61 . (18) 

To estimate the exponent r\ theoretically, we note that for 
predictability to be possible we need the divergent term 
of at least B2 to be of the order of unity. The quantity 
B2 (cf. Eq. (|12p) has a coefficient a which is of order of 
1/N because there is only a single localized mode that 
becomes unstable @, and the divergent term that goes 
like (7 — 7p)~ 3 / 2 - Thus we expect predictability when 
|7 — 7p| ~ iV~ 2 / 3 , which estimates 77 = —2/3. Higher 
order derivatives are more singular and thus r\ can only 



be smaller due to their effect. This estimate is very im- 
portant, since it guarantees that predictability will not 
deteriorate in the elasto-plastic steady-state. To see this, 
we recall that the average distance between successive 
plastic events (A7) is known [l2[ to follow a scaling law 
with the system size (A 7 ) ~ . The exponent /? was 
measured in a variety of systems in 2D, and was found 
to be always in the range [-0.67, -0.63]. In the present 
model system we measure (5 ~ —0.65. With the present 
accuracy we cannot exclude that 77 rs /3, leading to the 
realization that while the predictability is reduced with 
increasing N, the range of 7 over which we need to pre- 
dict reduces almost at the same speed, if not slightly 
faster. Thus in effect the predictability in the elasto- 
plastic steady state does not deteriorate with N. 

In summary, we showed that (i) the discussion of plas- 
ticity in amorphous solids calls for understanding the 
role of nonlinear elasticity, and (ii) for the elasto-plastic 
steady state the knowledge of a few nonlinear elastic con- 
stants suffices to predict where the next plastic event 
should occur. Finally, the analysis presented above is 
global, but it can be extended to local instability maps 
eliminating the system size dependence whatsoever. The 
first step toward a local analysis, which is the microscopic 
definition of the nonlinear elastic constants, has been al- 
ready achieved and is available in Ref. [l5j]. The second 
step where it will be shown how to use these results in 
models of elasto-plasticity in both two and three dimen- 
sions will be discussed elsewhere Hal. 
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